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Parameter Estimation for a partially observed Ornstein-Uhlenbeck process 

with long-memory noise 
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Abstract: We consider the parameter estimation problem for the Ornstein-Uhlenbeck pro¬ 
cess X driven by a fractional Ornstein-Uhlenbeck process V, i.e. the pair of processes defined 
by the non-Markovian continuous-time long-nremory dynamics dXt = —dXtdt + dVt ; t ^ 0, 
with dVt = —pVtdt + dB\ t ^ 0, where 9 > 0 and p > 0 are unknown parameters, and 
B h is a fractional Brownian motion of Hurst index H £ (^,1)- We study the strong con¬ 
sistency as well as the asymptotic normality of the joint least squares estimator (0t, Pt^J 
of the pair (#,p), based either on continuous or discrete observations of {X s ; s £ [0, T]} as 
the horizon T increases to +oo. Both cases qualify formally as partial-observation questions 
since V is unobserved. In the latter case, several discretization options are considered. Our 
proofs of asymptotic normality based on discrete data, rely on increasingly strict restrictions 
on the sampling frequency as one reduces the extent of sources of observation. The strategy 
for proving the asymptotic properties is to study the case of continuous-time observations 
using the Malliavin calculus, and then to exploit the fact that each discrete-data estimator 
can be considered as a perturbation of the continuous one in a mathematically precise way, 
despite the fact that the implementation of the discrete-time estimators is distant from the 
continuous estimator. In this sense, we contend that the continuous-time estimator cannot 
be implemented in practice in any na’ive way, and serves only as a mathematical tool in the 
study of the discrete-time estimators’ asymptotics. 

Key words: Least squares estimator; fractional Ornstein Uhlenbeck process; Multiple inte¬ 
gral; Malliavin calculus; Central limit theorem. 
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1 Introduction 

1.1 Context and background 

The question of drift parameter estimation for solutions of stochastic differential equations 
driven by fractional Brownian noise goes back at least as far as the seminar work of Kleptsyna 
and Le Breton in [T7] (also see prior references therein), where a maximum likelihood estimator 
(MLE) was proposed. This work was a first genuine attempt to show how to compute the 
MLE in practice in the regular case (self-similarity ‘Hurst’ parameter H > 1/2), by relying on 
continuous-time observation of a single path over a finite time interval, with strong consistency 
and asymptotic normality as the horizon increases to infinity. This work was itself motivated 
by - and is an extension of now classical ideas of how to use the Girsanov theorem to compute 
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the MLE in the case of white noise, as presented in the 1978 treatise [19] . The continuous¬ 
time data was also invoked in to justify that any diffusion-type parameter (any constant 
multiplicative term in front of the equation’s driving noise) would then be directly observable. 
That observation remains generally true in many contexts for continuous-time data, including 
when the noise is fractional Brownian. It explains why so many authors since ca have 
continued to study the estimation of drift parameters for fractional-noise-driven problems. 
We have listed some of these works below. Our paper inscribes itself in this line of work, and 
can be viewed as a study of partial observation questions, as we will explain shortly. 

Generally speaking, this effort to understand drift estimation for ergodic diffusions, even 
Gaussian ones, is of fundamental importance for any quantitative study where mean-reverting 
quantities are believed to be asymptotically stationary, and are either observed with noise or 
are intrinsically stochastic. The best known class of examples, which also encapsulates the 
question of whether or not the stochastic process of interest is observed directly or indirectly, 
is that of stochastic volatility in quantitative finance. Some of the original ideas on how to 
estimate this volatility’s drift parameters is given in the 2000 research monograph HZ]. Similar 
models with fractional noise were introduced in the late 1990’s, as in m for continuous time, 
but to our knowledge, their statistical estimation was left unexplored for more than ten years. 
The broadest question, applicable in finance and other fields, is to estimate jointly the drift, 
diffusion, and memory parameters for fractional-noise-driven equations; it is typically non¬ 
trivial, and with the exception of one study in [4], largely unresolved by bona fide statistical 
means for non-self-similar discretely observed continuous-time processes. We will not address 
this issue here. In the case of partially observed data, we refer the interested reader to a study 
in the context of high-frequency financial data, where a sequential Bayesian methodology is 
combined with classical estimation techniques and a calibration method to find H\ see [9]. 
We also mention a method described in m for a similar model with partial observation, 
where it is shown that the minimax-sense optimal estimator has a very slow convergence rate, 
and relies on high-frequency data; it can be argued that when applied to bona fide financial 
data, this estimator cannot be implemented without leaving the realm where continuous-time 
semi-martingale models with long-memory volatility are appropriate. 

The present paper looks specifically at a difficulty which arises when one single path is 
used to estimate more than one real drift-type parameter. Questions of identifiability can arise 
in this context (see for example the treatment of the Generalized Method of Moments (GMM) 
as described in EH). Once such a question has been resolved, the main practical obstacle to 
implementation is typically that of discretizing the data given by the path, i.e. using only 
observations in discrete time. For the sake of conciseness, we look at a specific situation 
where identifiability is resolved explicitly and in a natural way, within the task of deriving the 
continuous-time and discrete-time estimators, without having to rely on abstract conditions 
which would provide a priori identifiability of the vector of parameters. Our objective is to 
describe conditions under which strongly consistent and asymptotically normal estimation 
can be established quantitatively, based on discretely-observed data alone. Since access to 
and analysis of continuous-time data is not typically a realistic assumption, we will view any 
estimators based on continuous-time data as tools in the task of deriving strong consistency 
and asymptotic normality for the discrete-time estimators. Our strategy is to attempt to 
discretize the former, which includes a need to approximate Riemann and stochastic integrals. 
In the process, we discover some fundamental differences between the discrete and continuous 
time estimators. In particular, we find that while least squares (LS) estimation appears to 
be the best tool based on continuous-time observation of a path, when converted to discrete 
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time data, the estimator’s interpretation as a LS optimizer is lost, and a GMM interpretation 
seems more appropriate. This introduction, including our summary of results in Section Q221 
contains specific details supporting these ideas. 

We are largely motivated by the paper [6], which studies the drift-estimation problem 
for the Ornstein-Uhlenbeck process driven itself by another, unobserved Ornstein-Uhlenbeck 
process (OU-OU). Their work only deals with an underlying white noise. Specifically, let W 
be a standard Brownian motion and let 0 and p be positive real parameters. The OU-OU 
process is the solution of the following system 

f Xq = 0; dX t =—0X t dt + dVt, t ^ 0; , . 

1 U 0 = 0; dV t = -pV t dt + dW t , i > 0. U 


Here one may consider that X is observed, while V is not, or conversely, or that both processes 
are observed; which interpretation is used makes a crucial difference, as we will see. 

Since the quadratic variation of V is t, the classical Girsanov theorem implies that a natural 
candidate to estimate 9 is the MLE (recall that this idea was already contained in [19]), which 
can be easily computed for this model: one gets 


-fpXtSXt 

Jo X tdt 
-X$ + T 
T^xfdt 


( 2 ) 

( 3 ) 


In (]2]), the integral with respect to X must be understood in the Ito sense. Consequently, 
line ([3]) follows from Ito’s formula, and the fact that X too has quadratic variation equal to t. 
One also notes that this estimator’s construction is in fact non-dependent on the form of the 
bounded-variation part in the definition of V, which can be interpreted as a form of robustness 
of 9x with respect to model misspecification, although we are about to see that the behavior 
of 9x depends heavily on U’s drift specification. 

On the other hand, it is worth noticing that Ox coincides formally with a least squares 
estimator (LSE). Indeed, by interpreting j'J Xf Xf dt as the Ito integral in (]2]), Ox formally 
minimizes 



x t + ex t 


2 

dt. 


By reversing the roles of V and X , one can obtain an estimator for p similar to |3]). However, 
if V is unobserved, one may recast that estimator by using the estimated value of V) based on 
Ox and the observed path of X, based on dTJ. In other words, one defines 


—Vx + T 

PT — -- 

2 fo T V t 2 dt 

where V) = Xt + Ox Jq Xtdt for every t ^ T. As it turns out, this joint estimator (^0x- pr'j, 
based on continuous observations of X alone, does not converge to {0,p), but rather to an 
explicit rational function of the pair of unknown parameters ( 0,p ). This was proposed and 
proved in [6], wherein a semimartingale approach was used to study the asymptotic behavior. 
Specifically they showed 
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• strong consistency: as T 


+oo, almost surely, 


(4) 

(5) 


Ot 

Pt 


0 + p 

0p{8 + p) 

(6 + pY + 6p 


• asymtpotic normality : as T —> +oo, 


Vt 


(^t - {0 + p),pt 


0p(0 + p) \ 
(1 9 + pY + Op) 


law 


N(o,r) 


where T is a covariance matrix which has a explicit form as a function of 0 and p. While 
intuition gathered from the full-observation case is in fact erroneous when V is unobserved 
(the naive candidates for Ot and px lead to modified limits rather than (Oti'pt^J —> (0, p) ), 
nevertheless we have the full picture for the asymptotic behavior of the MLEs/LSEs associated 
with (HD, at the minor cost of having to solve a non-linear system of two equations to obtain 
consistent estimates of ( 6,p ). 

In the present paper, our goal is to investigate what happens when, in (jT]), the standard 
Brownian motion W is replaced by a fractional Brownian motion B H . Thus we assume from 
now on that X is an Ornstein-Uhlenbeck process driven by a fractional Ornstein-Uhlenbeck 
process V : this means the pair ( X , V) is the unique solution of the system of linear stochastic 
differential equations 


f Xq — 0; dXt — —OXtdt + dVt, t ^ 0 . . 

\ Vo = 0; dV t = -pV t dt + dBf, t ^ 0, 1 j 

where B H = {B^,t ^ 0} is a fractional Brownian motion (fBm) with Hurst index H G (4,1), 
and where 0 > 0 and p > 0 are unknown parameters. Though X is defined for all H in (0,1), 
to keep technical difficulties to a reasonable level, we restrict ourselves to the case H G (1/2,1). 
It turns out that we need the condition 0 Y P f° r identifiability; remarkably, this is not needed 
when H = 1/2, as we saw in the system (j4j), (|5|) . Details of our results are summarized in 
Section 11.21 


We now provide some references to estimation with fBm noise, which are further moti¬ 
vations for our work. We mentioned that the single-drift parameter estimation problem for 
fractional diffusion processes based on continuous-time observations was originally studied in 
m via maximum likelihood; more recent work on this question includes, e.g., [28( [26] , Re¬ 
cently, the LSE for the fractional Ornstein-Uhlenbeck (fOU) process, i.e. the process V in 
d6| was proposed in M- assuming V is fully observed in continuous time, the LSE for p is 
defined by 



where the integral J Q T VtdVt is interpreted in the sense of Skorokhod. This integral is the 
extension to fBm of Ito’s integral for Brownian motion. In the case p > 0, m proved that px 
is strongly consistent and asymptotically normal as T —> oo. Unlike in the case of Brownian 
motion (H = 1/2) discussed above, ~p T does not coincide with the MLE given in |17| . because 
the Girsanov theorem for fBm takes a different form than in the case H = 1/2. Given the 
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notorious fact that Skorohod integrals are difficult to interpret in practical terms for fBm, the 
authors of m proposed in addition the following alternate estimator, which is arguably a 
method of moments : 

G maw!,** ) (7) 

they proved it is strongly consistent and asymptotically normal. In the case p < 0, [5] 
established that Pt of p is strongly consistent and asymptotically Cauchy-distributed. 

The alternate choice of estimator (0) does not, however, avoid the use of continuous obser¬ 
vations over discrete ones; this is a problem with many works on fBm-driven models, and is an 
additional motivation for us to study the asymptotics of estimation for fBm-driven processes 
based on discrete observations. There exists a rich literature on this practical problem for 
ordinary diffusions driven by Brownian motions; we refer for instance to [26 j . A handful of 
authors are beginning to look at these questions with various fBm-driven models, starting with 
[28] in 2007, and more recently [3} 0 [U [20] • In particular, for the fOU process V, motivated 
by the estimator given in 0), [15] studied its natural Riemann-sunr-type time discretization 


( 1 

l nHT(2H) 



1 

2 H 


( 8 ) 


providing strongly consistency and Berry-Esseen-type theorems for it. While we have no 
doubt that this estimator is indeed strongly consistent and asymptotically normal, the proofs 
in [15] rely on a possibly flawed technique, since the passage from line -7 to -6 on page 434 
therein requires the condition H > 3/4, while one expects normal asymptotics only for the 
case H ^3/4. 


In our paper, we focus our discussion on estimators which are derived from a basic LSE, 
since that technique is known, at least in the Brownian case described in [6], to allow for a 
straightforward bivariate extension, as mentioned previously. The LSE has also given rise to a 
number of successful studies in the univariate case with fractional processes: we have already 
cited [3 ,0 TL [20], while [2j is the continuous-time version of [3]. 

Herein, specifically, we begin our study of LSE for ( 9 , p) in ([6]) by using the formal least- 


rj 

squares interpretation mentioned above, i.e. looking for the nrinimizer of 9 —> J Q 
this leads formally to the following estimator for 9 


x t + ex t 


dt 


Jo X t 6X t 


fo x ? dt 


(9) 


and to the similar estimator 



( 10 ) 


for p, where, because of the fact that V is unobserved, one uses Vt = Xt + 9 t Jq Xtdt for every 
t ^ T, instead of relying on the unobserved Vt in the construction of pr- These estimators 
9t, Pt are no longer the MLEs, since, as we mentioned, the Girsanov theorem for fBm does not 
have the same form as for Brownian motion, but they are still formally LSEs. Nevertheless, 
there is a major difference with respect to the Brownian motion case. Indeed, since the process 
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X is no longer a semimartingale, in Q and (fTOl) one cannot interpret the numerators using the 
Ito integral; the Skorohod integral turns out to be the correct notion to use here. We mentioned 
above that Skorohod integrals are difficult to use in practice, but since our Hurst parameter H 
exceeds 1/2, it is possible to reinterpret the Skorohod integrals as so-called Young integrals, a 
pathwise notion, modulo a correction term which we will be able to compute explicitly thanks 
to the Malliavin calculus. 

Having succeeded in correctly interpreting the stochastic integrals in (J9]) and (11011 , the issue 
of how to discretize them becomes paramount to practical implementation, and herein we will 
propose several different options, some of which allow for strong consistency and asymptotic 
normality under broader conditions than others. 

Our discrete-observation study also applies to the case H = 1/2 as a limiting case. The 
article [6j treats this case solely with continuous observations; our work thus covers an exten¬ 
sion of their work to discrete observations. Checking the validity of this statement rigorously 
is straightforward; for the sake of conciseness, we leave it to the interested reader. 

1.2 Summary of results and heuristics 

We now summarize our results, the structure of our article, and our main proof elements, 
including useful heuristics when available. 

• In Section 2 and in the Appendix we introduce the needed mathematical background 
material for our study, including elements of the Malliavin calculus, a convenient criterion 
for establishing normal convergence on Wiener chaos, and the relation between Skorohod 
integrals and Young integrals with respect to fBm when H > 1/2. 

• In Section 3, we concentrate on proving strong consistency and asymptotic normality 
for the estimators 8t and px with continuous observations. 

— We first prove the following almost surely convergences: 


8t 

—> 8* :— p + 8, 

(11) 

Pt 

* pd (8 + p) 

(12) 

" ' +(s+/>F 


The proof relies on studying the numerator and the denominator of the expres¬ 
sions for 9t and pr separately. For the denominators, we rely on Birkhoff’s ergodic 
theorem, and elementary covariance estimations for exponential convolutions with 
fBm. For the numerators, we express the Skorohod integrals as Young integrals 
plus their correction terms involving Malliavin derivatives which are explicit deter¬ 
ministic functions since our processes are Gaussian. 

— The expression 8* = p +8 in dill) is easy to explain: as noted below in line (1161) . X 
satisfies a stochastic integro-differential equation in which the zero-mean-reversion 
term — (p + 8) X (t) dt appears, and thus a natural candidate for a consistent esti¬ 
mator of p + 8 is the LSE 8, whether one adds merely one mean-zero noise term 
dB H or another term Xgds'j which is asymptotically small. This is why the 

limiting behavoir of 8 remains the same for us as in (H|) , which is the Brownian 
case (If = 1/2) studied in [6]. The details of this heuristic are omitted, since the 
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full proof we present herein is needed to be convincing. On the other hand, the 
expression in (112[) for our p* is more opaque; there does not seem to be a direct 
heuristic to explain it, beyond our computations. When one compares our p* in 
(fl2l) with the p* in (j5|) found in [6], one sees that the term pQ in Q is replaced by 
the expression (p 2 ~ 2H — 6 2 ~ 2H ^ / (0~ 2H — p ~ 2H ), which can help identify how the 
case of fBm deviates from the case H = 1/20 The expression for p* is analyzed 
further in the context of discretizing (j), , which helps explain to some extent why 

this complicated expression arises, as the reader will find out in the first paragraph 
of Section 14.21 


We prove asymptotic normality of ( 6t, PTJ by expressing the Skorohod integrals as 
iterated Wiener integrals, identifying dominant portions of these integrals, relying 
on a criterion for normal convergence in law in Wiener chaos, combined with a 
number of almost sure convergences. Our main asymptotics normality result is a 
central limit theorem that holds for \/T [Ot — 0*,pr — p*^j as T —>• oo, as soon as 
H € [1/2, 3/4). The asymptotic covariance is given explicitly. See Theorem[8l The 
upper limit of validity of this theorem is a typical threshold in normal convergence 
theorems in the second Wiener chaos. See for example a classical instance of this 
situation in the Breuer-Major central limit theorem, as presented in |22l Chapter 
7]. For H > 3/4, we conjecture that the estimators are asymptotically Rosenblatt- 
distributed (again see [221 Chapter 7] for a classical example of such a phenomenon), 
and that the convergence occurs almost surely; this point is not discussed further, 
for the sake of conciseness. 


The topic of Section 4 is to construct estimators based solely on discrete observations. 
The asymptotic results we prove still require increasing horizon. We also assume that X 
is observed at evenly spaced intervals, with a time step A n , and we set the time horizon 
to be T n = nA n —>• oo as n —>• oo. Let tk = kA n be the Lth observation time. For 
instance, the case A n = 1 corresponds to a fixed observation frequency; other conditions 
on A n will include requiring A n to tend to 0 as fast as a certain negative power of 
n, i.e. the observation frequency increases as the horizon increases. For some strong 
consistency results, it will even be possible for us to relax conditions on A n where it 
is allowed to tend to infinity like a power of n, i.e. with decreasing frequency as the 
horizon increases. 

Arguably, to be consistent with the assumption that only X is observed, the only estima¬ 
tors which are of practical use are those which rely solely on the values {Xt k : k = 1 ,,n}. 

Designing such an estimator by discretizing (Ot, Pt^J turns out to be a difficult task, in 
which the final expression solves a non-linear system in the spirit of that which would 
follow from dill) and (1 12[) , but is rather distinct from this system because of the difficulty 
in how to interpret the discretizations of the Skorohod or Young integrals. 

The method we have chosen moves through several intermediate steps, where we grad¬ 
ually increase the number of terms in the estimators which are replaced by discretized 
versions. This method has the advantage of clearly showing where the restrictions on the 


4 As a way to compare these terms, which do coincide when H = 1/2, we can see that if p tends to 6, 
the aforementioned expression in (j5j) tends to 9 2 (1 — H) /H, which thus deviates significantly from the case 
H = 1/2 quantitatively, particularly for H close to 1. 
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observation frequency A n come into play. Each intermediate estimator can be consid¬ 
ered as a perturbation of the previous one, starting with the continuous-time estimator 
of Section 3. Thus arguably all these estimators can be considered as tools used for 
the final objective, attained in Section 14.31 of constructing a strongly consistent and 
asymptotically normal estimator of ( 9 , p) based only on the data {X tk : k = 1,..., n}. 
Nonetheless, some of the other estimators are relevant in their own right, as they might 
correspond to realistic partial or full observation cases. 

— The main technical estimates which allow our discretization are given at the be¬ 
ginning of Section 4. These are Lemmas [TO] and [TT] based on applications of the 
Borel-Cantelli lemma. For Z any stochastic process, we let 


Q n (Z) ■= n~ l (Z tk ) 2 ■ 

k =1 

Let St ■= Jq Xgds and Ej := fg X s ds. We show that the discrepancy between 
ST n /T n and its discrete version Q n (X) is = o(l/\/T/) almost surely. We then 
compute three different discrepancies related to E: first we show that the difference 
between T" 1 J Q Tn E 2 cLt and its discrete version Q n (E) is also = o (1 /^/l/) almost 
surely. Then we show that with E the version of E which depends only on X 
observations, i.e. E t k '■= A n Y^i=i Xf i _ 1 , we get that Q n ^E^ — Q n (E) tend to 0 
almost surely. This is helpful to prove strong consistency of discrete estimators. To 
prove asymptotic normality, we need that Q n ^E^j — Q n (E) = o (l/\/7n")j which 
we prove holds almost surely. Increasingly restrictive conditions on A n are needed 
for these successive results. 

— We first concentrate on discretizing the denominators of (0x n , Pt, 

* We replace the denominator of 9t u by Q n (A), yielding an estimator 9 n , and 
we then replace the denominator of pr by Q n ( X ) + ( 9 n ) 2 Q n (E), yielding an 
estimator p n , because, as it turns out, V t 2 dt is asymptotically equivalent, 
almost surely, to Q n (X) + ( 9 n ) 2 Q n (E). Thanks to this, to the almost sure 
equivalence of Q n ( X ) with ST n /T n , and similarly for Q n (E), coming from Lem¬ 
mas m and EU the strong consistency and asymptotic normality of ( 9. 


n Pn 

follows from that of proved in Section 3. Here it is sufficient to 

assume H € (1/2,1) and A n ^ n a for some a € (—oo,l/iL) for the strong 
consistency; note that A n is allowed to remain constant or even increase like a 
moderate power in this case. For the asymptotic normality, it is sufficient that 
H € (1/2,3/4) and nA^ +1 -»• 0. 

* A second result is obtained in which we forego having access to the process E 
itself, relying instead on its discrete version T,t k = A n Yli=i X/._ ; in this case, 
almost-sure converge of p n requires n a+1 A ^ +1 — >• 0 for some a > 0, and the 
central-limit result for p n requires n 3 A„ H+3 —>• 0. 

We are then able to define and study a bonafide estimator based on discrete data 
alone. 


* We begin with assuming that we have access to both Xf k and Ej fc for all 
k = 1 ,n. The stochastic integrals in 0 n and p n were analyzed in Section 3, 
and were found, under scaling by T” 1 , to be asymptotically constant, where 
the explicit constants depend on the parameters. By using these limits and a 
discretization of the Riemann integrals in the denominators of 6 n and p n , this 
allows us, at the beginning of Section 14.21 to motivate the definition of a pair 
of estimators (0 n , p n ) as solution of the non-linear system 


F{O n ,Pn) = (Qn(X),Q n (S)) 

where F is a positive function of the variables (x,y) in (0,+oo) 2 defined by: 
for every (x, y) € (0, +oo) 2 


F(x,y) = HT{2H)x 


_^( y 2-2 H _ X 2-2H X -2H _ y -2 H) 

((1 — H)x~ 2H , F[x~ 2H ~ 2 ) if y = x, 


if x / y 


and the data statistics used in the system are Q n (X) and Q n (E). Strong 
consistency and asymptotic normality follow for the uniquely defined (0 n , Pn)- 
The delicate computation of the asymptotic covariance is given. The parameter 
restrictions remain the same as for (j ) n , p r ^j , namely strong consistency if H € 
(1/2,1) and A„ ^ n a for some a € (—00,1/#), and asymptotic normality if 
H € (1/2,3/4) and nA,f +1 ->• 0. 

* By redefining the estimators (6 n , p n ) using E instead of E, one ensures that 
only the data {Xt k : k = 1,.. , ,n} is used. Here, the results from the previ¬ 
ous case can be applied directly with the auxiliary results on how E perturbs 
E (Lemma HID , to obtain the same almost-sure convergence and central-limit 
result, but these are now restricted respectively to the aforementioned ob¬ 
servation frequency parameter ranges n“ +1 A^ +1 —> 0 for some a > 0, and 
n 3 A 2i?+3 —»• 0. 

— Finally, to illustrate how the complexity of the nonlinearities in the definition of 
(QniPn) may be attributable to the partial-observation problem, we define a pair 

of estimators (jt n ,p^j under the assumption that both {Xt k : k = 1,... ,n} and 
{Vt k : k = 1,..., n} are available. The p^ is explicit given {Vt k : k = 1,..., n}, and 
is identical to the one given in mi i.e. ®, as it should be. The 9_ n satisfies the 
following straightforward non-linear equation given p^and {Xt k : k = 1,,n}: 



The parameter restrictions are the same as when X and E are discretely observed: 
strong consistency holds for (j[ n , p^j if H G (1/2,1) and A n ^ n a for some a € 
(—oo, 1 /H), and asymptotic normality holds if FI € (1/2,3/4) and nA^ +1 —>■ 0. 


Before we proceed with the details of our study, we provide needed mathematical tools in 
Section 2. 
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2 Preliminaries 


In this section we describe some basic facts on the stochastic calculus with respect to a 
fractional Brownian motion. For a more complete presentation on the subject, see |23| and 

HI- 

The fractional Brownian motion (Bf 1 ,t ^ 0) with Hurst parameter H E (0,1), is defined as a 
centered Gaussian process starting from zero with covariance 

R H (t, s ) = E(B?B^) = i (t 2H + S 2H - 1 1- s \ 2H ) ; s, t ^ 0, 

We assume that B H is dehned on a complete probability space (H, T, P) such that T is the 
sigma-held generated by B H . By Kolmogorov’s continuity criterion and the fact E ( Bp 1 — Bf\ 
|s — t\ 2H , we deduce that B H admits a version which has Holder continuous paths of any order 

7 <H. 

Fix a time interval [0,T]. We denote by 'H the canonical Hilbert space associated to the 
fractional Brownian motion B H ; the book [23] . among many other references, can be consulted 
for the construction and properties of %. We use the following convenient notation for Wiener 
integrals with respect to B H : 

B h {ip) := I p(s)dB H . 

Jo 

Of interest to us is the fact that, with H > 1/2, for a pair of (non-random) functional elements 
ip, ip of "H, its inner product satisfies 

(ip, ip) n = E [B h (ip) B h (ip)) = H{2H — 1) f f (p{u)ip(v)\u — v\ 2H ~ 2 dudv. 

Jo Jo 

It follows from [25] that the set \H\ of functional elements in is Banach and actually contains 

L*([0,T\)- 

The Malliavin derivative D w.r.t. B H , which is an "H-values operator, is dehned hrst by 
setting that 

DB h (<p) = <p 

for any ip € H, and then by requiring that it satisfy a multi-parameter chain rule: for any / E 
C£°( M n , R) (inhnitely differentiable functions from M n to M with bounded partial derivatives) 
and any tpi, ...,tp n E "H, D operates on the cylinder r.v. F := f(B H (ip\),B H (p n )) as 


n 


DF = V— (B 

“ dxi 
1=1 


H 


{iPl ),-,B H (ip n ))iPi- 


The domain D 1,2 of D is then the the closure of the set of cylinder r.v.’s F with respect to 
the norm 

Ill'll 1,2 :=E(F 2 ) + E(\\DF\\ 2 n ). 

The divergence operator 5 is the adjoint of the derivative operator D : an ^-valued r.v. 
u E L 2 (H;"H) belongs to its domain DomS if 


E\(DF,u) n \ ^ c u ||F|| L 2 (n ) 


10 


for some constant c u and every cylinder r.v. F. In this case 5(u) is uniquely defined by the 
duality 

E(F5(u)) =E(DF,u) n 

for any F £ D 1,2 . We will make use of the notation 



u £ Domd. 


In particular, 5 extends the Wiener integral: for h £ \TL\, B H (h) = 5(h) = | Q T h s 5B^. 

For every n ^ 1, let 'H n be the nth Wiener chaos of B H , that is, the closed linear subspace 
of L 2 (Ll) generated by the random variables {H n (B H (h)),h £ TL, \\h\\u = 1} where F[ n is the 
nth Hermite polynomial. The mapping I n (h® n ) = n\H n (B H (h)) provides a linear isometry 
between the symmetric tensor product 'H® Tl (equipped with the modified norm ||.||-^©n = 
^tII-IIh®") an d Hn- It also turns out that I n (h® n ) is the multiple Wiener integral of h® n 
w.r.t. B H . For every f,g £ E &n the following product formula holds 


E (In(f)In(g)) =n ] .(f 1 g) H ® n . 

For h £ E® n , the multiple Wiener integrals I q (f ), which exhaust the set E q , satisfy a hyper- 
contractivity property (equivalence in E q of all L p norms for all p ^ 2), which implies that 
for any F £ we have 

(E[\F\P]) 1/P ^c m (E[\F\ 2 ]) 1/2 for any 2. (13) 

It is well-known that L 2 (Q) can be decomposed into the infinite orthogonal sum of the spaces 
E n . That is, any square integrable random variable F £ L 2 (Q) admits the following “Wiener 
chaos” expansion 

OO 

F = E(F) + J2Ufn), 

n =1 

where the f n £ TF' )n are uniquely determined by F. 

Finally, we will use the following central limit theorem for multiple stochastic integrals (see 

PH). 


Theorem 1 Let {F n ,n ^ 1} be a sequence of random variables in the q-th Wiener chaos E q , 
q ^ 2, such that lim^^oo E(F 2 ) = a 2 . Then the following conditions are equivalent: 

(i) F n converges in law to Af(0,a 2 ) as n tends to infinity. 

(ii) converges in L 2 to a constant as n tends to infinity. 


3 Asymptotic behavior of LSEs 

Throughout the paper we assume that H £ (|, 1), 9 > 0 and p > 0 such that Of- p. 

It is readily checked that we have the following explicit expression for Xp. 

= (14) 


11 




where for m > 0 


xr = 


f 


0 -m(t- S ) dB H_ 


(15) 


On the other hand, we can also write that the system ({6]) implies that X solves the following 
stochastic integro-differential equation 


dX t = -(8+ p) X t dt - p9 Qf X s ds^J dt + dB?. 


(16) 


For convenience, and because it will play an important role in the forthcoming computations, 
we introduce the following processes related to Xp. 

S T = f X 2 dt ; £ t = f X t dt ; L T = [ V 2 dt~ P T = [ X t V t dt ; 

Jo Jo Jo Jo 

f 


and 


Lt = 


V 2 dt 


where for 0 ^ f ^ T 

V t = X t + 0 r E t , (17) 

and 8t is our continuous LSE for 0 as given in ((9|). We will need the following lemmas. 

Lemma 2 Assume H E (^, l). Then, as T 


oo 


X^dt —> r / x , 


almost surely, where 


and 


- t j: t 

j; fo E tdt — 

1 f T 

- / Z f X f dt —> 0 

1 Jo 


r, x = ™ \p 2 - 2H - 8 2 ~ 2H ] 


(18) 

(19) 

( 20 ) 


p 2 -8 2 


= H l (2H J [8~ 2H ~- 2H i 


p 2 — 8 2 




Proof. From ([6]) we can write 

d 


X t 


= A 


where A = 


1 + p -dp 
1 0 


X t 
St 

'X 


dt + d 


B? 

0 


. The process geometrically ergodic because the largest 

eigenvalue of A is negative. Then to prove Lemma (2j using Birkhoff’s ergodic theorem (for 
instance see mi it is sufficient to study the convergence of E[X 2 ], E[E 2 ] and E[Et_Xt] as 


12 




t —> oo. 

For the convergence of E[X 2 ], (1141) leads to 




(^) 2e [W) 2 ] + (^) 2 e k^) 2 ] 


2 Op 

¥~W 


E[(x e ) t (xn}. 


Since 

vX = {j^d) x{p ’ l,)+ {i 37 ) 

then by using 1) of Lemma [21] we obtain 


2 Op 


A (M 


|r/ x -E[X t 2 ]| <c(iL,0,p) e -*/ 2 . 


Thus we deduce the convergence (1181) . 
Using the same argument and the fact that 


St = 


0 


X? - xf 

p-e 


( 21 ) 


we deduce the convergence (fl9l) . 

Finally, the convergence (1201) is satisfied by using f 0 T^X^dt = ^ and point 5) of Lemma I2T1 


Lemma 3 We have 


± £ XtdXt—y-ip + d)^ 


almost surely as T —> 00 . 

Proof. From (jS]) and (IA-21) we can write 

r T r T rT 


f X,dX, = -e f Xfdt - p f XtVtdt + [ X t 5Bf 
Jo Jo Jo Jo 

= -0 [ X 2 dt -Op [ XtE t dt - p [ X 2 dt + [ X t dB{ 

Jo Jo Jo Jo 

r T r t 


an f f D s X t {t — s) 2H 2 dsdt 

Jo Jo 


where an = 2H(2H — 1). Moreover, 


f 

Jo 


X t dB, H = 


X 2 


X t dX t + (0 + p ) [ X 2 dt + Op f X t T, t dt 

Jo Jo 

+ (0 + p ) [ X 2 dt + Op [ X t E t dt. 

Jo Jo 


Thus 


pT v2 pT pt 

/ X t SX t = ^-a H / D s X t (t - s ) 

Jo 2 Jo Jo 


\2H-2 


dsdt. 


( 22 ) 
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Since 


D s Xr = e-^- s h m {s) 

we deduce that 

k.[ Xt6Xt = jf ^^(pe-^-ee-^ii-sfn-^dsdt 

- <23) 

Thanks to l’Hopital’s rule, as T —> oo 

^ [ T f J_( pe pr_ de 0r\ r 2H-2 drdt ) 2 _ 2H _ q2 _ 2H 

T Jo Jo P ~ 9 V ) p-6 

(p + 0)rj x 
olh 

Finally, combining this last convergence and point 5) of Lemma I2T1 the proof of Lemma [3] is 
done. ■ 

We now have all the elements to obtain our strong consistency result for 9t- 
Theorem 4 We have 

e T —> e* 

almost surely as T —> oo, where 9* = 9 + p. 

Proof. The proof follows directly from the convergence (1181) and Lemma 02 ■ 

The next lemmas are additional elements needed to prove the strong consistency of px- 

Lemma 5 We have ^ 

h: -v + ( „ + «) v 

almost surely as T —> oo. 

Proof. The equation (1171) ensures 

pT pT pT 

L t = X?dt + 29 T / X t Titdt + 0|i / T? t dt, 

Jo Jo Jo 

and the desired conclusion follows by using Lemma [2] and Theorem [4j ■ 

Lemma 6 We have 

j, J VtSVt —* -p9(p + 9)if 

almost surely as T —> oo. 

Proof. From (j6]) and (fTfl) we can write 

[ v t 5v t = [ x t sx t + e T [ x?dt-e T (9 + P ) [ z t x t dt- P 99 T [ s t 2 dt 

Jo Jo Jo Jo Jo 

f T r T 

+ / T, t 9 T 6B^ + / X t X t dt. 

Jo Jo 
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On the other hand 
f T 


/ 
J o 


Or^tdB, 


H 


f dr'EtdVt + p0T [ 
Jo Jo 


= -e 

= -9 


-7 

\L 


0 T X t dV t + 9 

T r T 


l 0 T [ V t dV t + p9 T [ Ti t X t dt + P 99 T [ Y? t dt 

Jo Jo Jo 


9 T X t dX t -9 T / Xfdt + 0- 
/ 0 Jo 


X t dt + ^ V T + P 9 t j T, t X t dt + p69 T j Y? t dt. 


[ X 2 dt + 9~ l 9 T [ 

Jo Jo 

■f 


V t dV t + P 9 t / ZtXtdt + P 99 T / Yi^dt 

o Jo Jo 


7 


= 2 7^/ 

Now, applying (1A-2|) . we obtain 

£ VtdVt = £ X t 5X t - 7_x! + e _Ty2 _ q t q £ z t X tdt + 9 2 t £ E t X t dt 


~otH [ [ D s (9T^t){t — s) 2H 2 dsdt. 

Jo Jo 


(24) 


On the other hand 
It follows from (1231) that 
9t = 

Hence, for s < T 


D s (9 T Y,t) = Y, t D s 9 T + 9 T D s Y, t . 


Jo X t S Xt 


S T 

\ x t ~ Jo Jo ( P e ~ pr ~ 9e ~ 6r ) r 2H ~ 2 drdt 

S T 

XtD s Xt — 9tD s St 


D s 9 t = 


St 


Thus 


,-T rt _ 

an D s (9T'Sit){t — s) 2H ~ 2 dsdt 



o Jo 


X T 


T rt 


9t 


T rt 


= I I X t D s X T (t-s) dsdt — a H Q 

JT Jo Jo ^T 

+o.h9t [ [ D s T, t (t — s) 2H ~ 2 dsdt 

Jo Jo 

:= JgT — -h ; T + Jo,T 


T, t D s ST{t — s) 2H 2 dsdt 


o Jo 


We shall prove that for every e > 0 

\Ji,t\ _ 

~e 

\ J 2,t\ 


0, 


0, 


T £ 

Jq,t _ 9 + p 

T p — 9 


HT{2H)[(- P y- 2H - (- 9y- 2H ] 


(25) 

(26) 
(27) 
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almost surely as T —>• oo. 

We first estimate Jit- Clearly, (1141) implies 


J\,T = 


a n-^~ [ [ T t D s X T {t - s) 2H 2 dsdt 

£>t Jo Jo 

^ (pe- p{T ~ s) - 0e" e(T " s) ) (t - s) 2H ~ 2 dsdt 

J ^tj (pe- p{ ?- t+x ^ - 6e - e{ ?- t+x )) x 2H ~ 2 dxdt. 


oh Xt 
p — 6 St 
otH Xt 


p — 0 St 

The last equality comes from making the change of variable x = t — s. 

\X t \/T £ sup te[0)T] || 


Hence 


\Ji,l 


c(H,9,p)- 


T e S t / T T e 

Using (1 18 [) . (1211) and the point 5) of Lemma |2H the convergence (l25l) is obtained. 
Next we estimate J^t- By (El) we have 


'h.T = OLTL 


9t 

!fr 


T rt 


II 


X t D s S T (t - s) 


2H-2 


dsdt 


2a H ^~ [ [ St [ X u D s X u (t - s) 2H 2 dudsdt 

jt Jo Jo Js 


= 2a H 


Ot 


T rt 


II 


Then 

\J2,t\ 

T £ 


< c(H,e,p) 

^ c(H, 9,p) 
—> 0 


(p - 9 )S t Jo Jo 
9t 


X„ 


J, J, J ' u 


(i - s) 


2H-2 


dudsdt. 


su Pte[o,T] \ x t\ su Pte[o,T] l^tl 


[ T f e ->"in(M(T-')(i _ s ) 2H ~ 2 dsdt 

Jo Jo 


S T 

9t 


T e 


su Pte[0,T] l^f|SUPt e[0jT ] |2Jt| 


St/T 


T e 


almost surely as T —> oo. The last convergence comes from (1181) . (I21j) . Theorem 0] and the 
point 5) of Lemma I2T1 Thus, the convergence (1261) is satisfied. 

Finally, we estimate Jst- Using (1211) and m 


q.h9t [ [ D s ll t (t — s) 2H 2 dsdt 

Jo Jo 


J3,T = 


r * 
P~ 8 J o Jo 


aH0T Jo Jo ( e ~ 6(t ~ S) ~ e ~ P[t ~ S) ) ~ s f H ~ 2dsdt 


By l’Hopital rule we obtain 


■h: 


9 + p 

~P~~8 


HT(2H){9 


1-2H _ pl-2H j 


almost surely as T —> oo. 

Using the above estimations fj 25 [) . (126[) . (127|) together with (1241) . Lemma [2J Theorem [4] the 
point 5) of Lemma [2D and Lemma [3] the desired result is then obtained. ■ 
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Theorem 7 We have the almost sure convergence 

Pt —> P* 

as T —>• oo, where 

= 9p(9 + p>y s 

r] x + (9 + p) 2 ry s ' 

Proof. The proof is a straightforward consequence of Lemma [5] and Lemma [6J ■ 

Our approach to prove the asymptotic normality for both estimators 9t and pt looks first 
at the normal convergence of the T-indexed second-chaos sequence based on the kernel which 
appears in the representation (1151) of X. Thereafter, thanks to elementary stochastic calculus 
in the second chaos, these double stochastic integrals will be identified in an expression for 
the leading terms in 9t — 9* in the proof of Theorem [8j A similar technique, plus the use of 
the chain rule of Young integrals and their relation to Skorohod integrals, is used to find again 
that the leading terms in pt — p* are also linear combinations of the same double integrals; 
the analysis of the lower-order terms are less evident than for 9t — 9*; the proof of Theorem 
[8] records all the details. 

We are ready to prove the asymptotic normality of ( 9t,Pt ) • 


Theorem 8 Assume that H € (|, |). Then 


Vf (e T - e*,p T - p*) ^ AA(0, 4 PT P) 

where the matrices T and P are defined respectively in L4-5I) and (,9,511 . 

Proof. We express 9t — 9* and pr — p* as linear combinations of the double stochastic 
integrals identified in the previous theorem, plus lower-order terms. 

The case of 9t — 9* is rather straightforward. It follows from (1161) that 

g T _„, = 

St 


Since 


and 


1 

vf 




XfYjfdt — — X 


2 

T 


1 

(, P-9)VT 

1 

(p-9)VT 


[ [ (pe"^" s) - 9e- e d~V5B^5B? 

Jo Jo 

(pJ 2 (/£)-0/ 2 (/£)) 


we can write 


where 


Vr(e T 


-9*) = 


w {eW e T )-Ph{m 


St/T 


+ R 9 t 


, e = P 9 [VT 

1 ' 2 S t /T 


(28) 
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(29) 











almost surely as T —> oo. 

For 'pi — p *, the situation is significantly more complex. We have for every 0 ^ t ^ T 

v t = x t + e T E t = v t + (0 T -e)E t = Vt + (e T -0*)x t + pE t 
= Vt-- 9 (X t -Vt)-^(e T - e*)(x t - v t ) 

= e l Vt -^x t - l -(e T -e*){x t -v t ) 


which leads to 

r T ^ 

Lt = / = It + (Qt ~ 9*)(Jt + (9t — 9*)Kt) 

Jo 

where 

It = -^{p 2 St + {0*) 2 L T - 29*pPr), 

Jt = -^{2pST + 26* L t — 2(9 + 2p)P T ), 

Kt = -^{St + L t — 2P t ). 

Thus, 

L t {pt ~ p*) = It+ ( 9t - 9*)(J% + ( 9 t - 9*)K%) 

where 


T v 

i T 

4 


f 


VtSVt - P*I T 


~P*Jt, and K^f = 


V 2 f T r t — 

= -^- + <* H / D s V t (t 

2 Jo Jo 

-P*K t . 


s) 2H 2 dsdt — p*It , 


On the other hand, using the formula (1A-1D we obtain 


’ S T 

< Pj, 

Lt 


iLz _i_ 
20 * ' r 
XtVt 


e * jo 
+1 


fo X s dB«- 


- f T V dB 
2p + p Jo v sun 


_p0_y2 
20* ^Ti 

f T y jdA i Yx 

Jo yi s uu s t 20* J 


H 


p JO 


Furthermore, using the relation between Young and Skorohod integrals, 


r S T 

< 

Lt 

< 


+ W fo X - SB . H + 7- Jo T Jo D.Mt ~ 
+ X,6B» + ff / 0 T f‘ D,X,(t 

- V 4 + 1 fl V - SB ? + f g So D >U(f - ») 


s)2H~ 2 d S dt — 

— s) 2H ~ 2 dsdt + 
2H ~ 2 dsdt. 


Setting 


A t ’■= OiH 



e-J- s \t-s) 2H ~ 2 dsdt 


(30) 
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we can write 


A* := 


:= 

A T ■= 


Ay := 


a H [ [ D s X t {t — s) 2H 2 dsdt = (p 2 2H - 6 2 2H ) X T 

Jo Jo p — 9 x 


T r t 


a H 


D s Vt{t — s) 2H 2 dsdt = p l 2tl Xt 


1 - 2 H ' 


1 0 Jo 

rT r t 



a H I I D s Yi t (t - s) 2H ~ 2 dsdt = \{Xl - Ay) = ( 6^ - p 1 - 211 ) X T 

'o Jo u p — v 

otn [ f D s Vt(t — s) 2H 2 dsdt = Ay + 0yAy + — J 2 ,t 

Jo Jo 

= Ay + 6* Ay + (0T — 0*)Ay + Jl.T ~ J2,T- 


The last equality comes from the fact that D s Vt = D s Xt + OxDsYit + SjZ) s 0y. 
Since 

\X , n* \T, _ P* ( P 2 \X (^*) 2 

x T+ e x T - ¥ y--x T ~ — 

we can write 


Ay + 2pAy 




V 2 

I¥ = -f + %-p'It 

= (e T -d*) a?-^ 


„2 rT /a*\2 rT 

MP+y)/ *.«*? +4 

u JO P JO 


j 

JO 


V S 5B S 


H 


+ -Rr 


where 


—V 2 

R T = Vt 


2 


21 X T P°r 2 , fl *x2^ on*„,- X T V T , 4 1 




H-R) 

2(9* 7 


Combining previous estimations we obtain 


Ry 


Lt(pt — p*) — CyRl/y) + Cy/2(/y) + —(^Ay + Jy + (0y — 9*)K T J + Ry 


where 


Cy = 


—pAy pp*(-2p + fc) p*(0*) 2 pjv p(0 T -0*)KJf 


(, p-6)S T 9 2 (p — 9) p6 2 ( p-6)S t ( p-6)S T 

-pA s pp*(-2p+^) p*(r ) 2 pA J 


= 


(, p-e)p x 9 2 (p — 9) 

almost surely as T —> oo, and 


pO 2 (p — 9)rj x 


(31) 


4 = 


9X £ 


| ^*(- 2j0 + ^) p*(r ) 2 | 


0jJf 


(p-0)Sr 9 2 {p-9) 


+ 


9(9 t - 6*)RV 


p9 2 ( p-9)S T (p-O)St 


c 9 = 


9 A- 


(p - e W 


| g p *(-2p + ^) p*(r) 2 | 


0A J 


P 2 (p-0) 


p$ 2 (p — 9)r] x 


(32) 
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almost surely as T 


oo. These last two convergences come from the fact that 


and 


A|/T —► A e = 2ffr(2 5 } - p 1 - 2 ^] 

p — 6 


Ji/T —> A ' 1 = -A [0*P - HT(2H)9*p- 2H ] 


almost surely as T — > oo, because A x/T — > 2HT(2H ) as T —>• oo. 
Thus 

7 F {c P tWt) + C tWt)) p 

+ jttrp 


VT(pt-p*) = 


~T~ 


where as T —> oo 


R? = R e T 


At + Jt + J i? T /VT 


Lx 


Lt 

T 


(33) 


(34) 


almost surely. 

Finally, with the expressions (j28|) and (13311 on hand, and the almost-sure negligibility of their 
corresponding lower-order terms as proved in (129[) and (1341) . we get 


Vr (o T - 6\px - p*) 



/ e t 
I p-e St 
P T 

\e-pS T 



+ (Rti K) 


where as T 


oo 


/ e t 
/ p-e St 
P T 

\e-pS T 



P : = 


e i 

p—9 Tj X 
P 1 

6 — p T} X 



(35) 


almost surely. Now, applying Slutsky’s lemma and Theorem [22] combined with the above 
convergences, the proof is complete. ■ 


4 Discrete observation 

Assume that the process X is observed equidistantly in time with the step size A n : ti = 
iA n ,i = 0 ,n, and T n = nA n denotes the length of the ‘observation window’. The goal 
of this section is to construct two estimators 9 n and p n of 6 and p respectively based on the 
sampling data Xf. , i = 0 ,n, and study their strong consistency and asymptotic normality. 
We also want to define estimators in such a way that consistency and normality results proved 
in Section 3 for the continuous-data estimators (Ot, Pt^J can be used to good effect in the 
discrete case. The basic strategy for this is therefore to look for ways of discretizing the MLE 
studied in Section 3. It turns out that the most efficient way of implementing this strategy 
is to define several intermediate estimators, starting with ones where only the denominators 
in ^ Ox , px^j are discretized, and then using an algebraic asymptotic interpretation of the 
numerators to avoid a direct discretization of the corresponding Young or Skorohod integrals. 
This method allows a rather direct use of the asymptotic normality Theorem [8] in Section 
3, while for the strong consistency results, some of the almost-sure convergences proved in 
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Section 3 are used directly, and additional ones are newly established early on in this section. 
See Section P for other details about the heuristics which explain the choices made below in 
this Section. 

For any given process Z, define 

1 n 

Qn( z ) :=-E ( Z ^) 2 

i= 1 

The following well-known direct consequence of the Borel-Cantelli Lemma (see e.g. CHI), will 
allows us to turn convergence rates in the p-th mean into pathwise convergence rates. This is 
particularly efficient when working with sequences in Wiener chaos. 


Lemma 9 Let 7 > 0 and po E N. Moreover let ( Z n ) n be a sequence of random variables. 
If for every p ^ po there exists a constant c p > 0 such that for all n E N, 

(E\Z n n 1/p ^c p -n-\ 

then for all s > 0 there exists a random variable r/ £ such that 

| Z n \ ^ r] £ • n -7+£ almost surely 


for all n E N. Moreover, E\r] £ \ p < 00 for all p ^ 1. 


As before we assume that A n = tk +1 — tk is a function of n only. Of some importance, 
particularly for the purpose of proving normal convergence theorems, is the case n~ a with a 
given a E M. The case a > 0 implies that the observation frequency must increase even as the 
horizon itself also increases. The case of a = 0 is of special importance because it corresponds 
to a setup where the observation frequency is fixed ( A n = 1, no in-fill asymptotics, only 
increasing horizon), which may be desirable in some applications. We will see that for some 
almost-sure convergence results, we may even take a time step A n which grows with n. In 
other words, this allows for very sparse observations. We will also see that most almost-sure 
results are valid for the entire range H E (^, 1), while normal convergence results require 
H E (^, |). We begin by recording and proving some important technical estimates. 


Lemma 10 Define 5 n (X) := s/Tf ^ Q n (X ) — ■ Then 

E [<£(*)] < c(H , e , p) min (nA 2 n H+ \ -i- + A ^ +1 + A ^' 3 £ f H ~ A j 

In particular, if A n ^ n a for some a E (— 00 ,1 /H), then 

Qn(X) —► 

almost surely as n —>• 00 . 

Proof. The points 3) and 6 ) of Lemma [21] lead to 

n fti rtj 


e K(v)] = te r r e 

n i j =1 ^ 1 J tj —1 

< c(H,e,p)nAl H+1 . 


x 2 t-x 2 \ (xl-x 2 ^ 


dsdt 


(36) 


(37) 
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On the other hand, we have 

E [<£(*)] 

-i n rt nt 

w ^ / L% l L o 



E (x 2 - X 2 ^ ) E (X 2 S - Xl , ) dsdt 


3 -1 


n i,j=l 1 ^ti-i 

i n rU rtj 

+Y E / / E [( x * - ^-i) - x *i-i)] E t( A * + x u-i) + x h-i)] “ 


n jj=i J U-\ Jtj -1 

n rU r l i 


+ 



^ , E [{X t - X ti _J (X, + X^.J] E [(X t + X^) (X s - X^)] dadi 

J -n 

1 

rE (■ D l(*>j) +- D 2(*,j) +£>3(*,j)) ■ 


L n • ■ -i 

*,j=i 


By using the points 2), 4) and 6 ) of Lemma [21] we obtain 


n 

— Y Di(i,j) 


i,j =1 


1 

T n 


n ntn 

£/ 

,i=i • y *»-i 


2 


E [ X? - Xl , ) dt 




€ 


c{H,e,p) 

Tn 

c{H,e,p ) 


-i 2 


A «Zl e ti_l/2 


i =1 

A„ 


^ g— An /2 


1 71 


i 


2— 1 

71 i /-t* 


, E., , 

n i = 1 ti— 1 J ti— 1 

< c(F, 0 ,p)A^ +1 , 



2(£(X*X S )) 2 + 2(£(X 2 _ 1 )) 2 - 4( J L(Xi i _ 1 X s )) 


dsdt 


and 


1 n 

- ]T (Z) 2 (i,j) + £>3(i,j)) 


*/i=i 




n *#j=i 

c{H,e,p) 


i rU 



U -1 ^ ti — l 


[( E(X,X,)f - (E(X u _,X,)f - (E(X ti _ 1 X t )) 2 + (E(X t ,_,X tj _,)) 2 ] dsdt 


E A‘"- 2 |i-i-l| 


4iT-4 


*<j=l 


4H-4 


< c(iL,0,p)A^- 3 ^j 

i=l 

Thus (1361) is obtained. Now, using (1131) . Lemma H and (1181) we will be able to assert the 
convergence (EZD, and thus the entire lemma, as soon as we can show that the right-hand side 


22 














of (1361) divided by T n converges to 0 as fast as some negative power of n. Thus we only need 
to show that there exists e > 0, such that as n-> oo 


Qn ■= mm 


A 


2 H 


(nA n )' 


+ ^l + -A 4H ~ 4 
n n 


Y.J 

3= 1 


4H-4 


^ n 


Let us concentrate first on the second part of the minimum defining q n . This is the sum of 
the three ter m s (nA n ) -2 , A %/n, and n^A 411-4 ^" =1 j 4H ~ 4 . The first of these three terms 
will tend to 0 like a negative power of n as soon as there exists £\ > 0 such that A n ^ n~ 1+ei . 
The second term will tend to 0 like a negative power of n as soon as there exists £2 > 0 such 
that A n ^ n 4 / H ~ £2 . For the third term, we must separate the case H < 3/4 from the case 
H ^ 3/4. When H < 3/4, the series J/” = i j 4H ~ 4 is bounded, so the last term in the second 
part of the min in q n will tend to 0 like a negative power of n as soon as there exists £3 > 0 
such that A n ^ n _1 /( 4_4 ^)+ £ 3_ When H > 3/4, the series is bounded above by a constant 
times n 4H ~ 3 , yielding a contribution of (A n /n) 4H ~ 4 ', so the last term in the second part of 
the min in q n will tend to 0 like a negative power of n as soon as there exists £4 > 0 such that 
A n ^ n -1+£4 . The case H = 3/4 is done in the same fashion, with the same conclusion as 
when H > 3/4. Thus we have proved that for each fixed n, if there exist £ 4 , £ 2 , £3 > 0 such 
that 

max (n _1+£1 ,n~™ +£3 ) ^ A n ^ n 1/H ~ £2 (38) 

then for some £ > 0, 

q n < n~ £ . 


On the other hand notice that for every H G (1/2,1), there exist £1,63 > 0 such that 


max (n 1+£l ,n 4-4/r+ £ 3^ '/ 2 _ (39) 

Thus for each fixed n, if we have 

A n ^ n~ 1/2 , (40) 

using the first part of the min in the definition of q n , we get 


q n ^ n 


with £ = H . To conclude, by (139[) . for each fixed n, we are either in the case (f38l) or (140p . so 
that q n ^ n~ £ in all cases as soon as A n ^ n 4 ' H ~ £2 for some £2 > 0. The proof of the lemma 
is complete. ■ 


Lemma 11 


Define <5 n (E) := \fTfi ^Q n (S) — f^ n E 2 dt^j . Then 


E [4(E)] ^ c(H, 9 , p) min [ nAf+ J , -±- + A^ +1 + A 4 "' 3 £ 3 

71 j =1 


4H-4 


In particular, if A n ^ n a for some a G (—00,1 /H), then 

Q n (E) —> if 


(41) 


(42) 
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almost surely as n —>• oo. 

On the other hand if n l+a A^ +1 —>• 0 for some a > 0, 


almost surely as n —>• oo, where 


Qn(S)-Q„(S) 


i=i 


(43) 


In addition, if n 3 A 2 2 H+3 —> 0, 


m L 2 {0) as n —>• oo. 


\/ik Qn(S)-Q n (S) 


(44) 


Proof. By using same arguments as in the proof of Lemma HOl (12111 and (1191) . we obtain 
(1411) and (l42j) . 

Now, we prove the convergence (1431) . We can write 


-1 


1=1 


Q n (E)-Q n (S) = — 2 

2—1 


2 2 


Using the point 6) of Lemma [2T] 


£ I ( S, , - E* 


EE/ / 

3=1 fc=l "'L-l Jt k-1 


-X tj _i)(X r -X tfc _i)]drd S 


I 1 1 /'L 

< c(i/, 0 ,p) x; / 

V 3=1 Jt ^ 1 


1 5 — ds 


< c(H,e,p) (nAf +1 ) 2 

Then, by Holder inequality and the point 3) of Lemma I2T1 we obtain for every p ^ 1 

PI \ 1 /p 

£ [ Qn( E) - Q„(S) J ) c(fd, 0, p) [n 2 A 2H+2 + nA* +1 ] . 

Thus, by (1131) . Lemma [9] and that fact that n 1+a A(/ +1 —>• 0 for some a > 0 the convergence 
(1431) is obtained. 

Furthermore, it is also easy to see that the convergence (j44h is satisfied. ■ 

4.1 Auxiliary estimators 6 and p 

The first step in constructing a discrete-observation-based estimator for which the asymptotics 
of (QtiPt^J studied in Section 3 can be helpful, is to consider the following two auxiliary 

estimators 9 n and p n of 9* and p* respectively, by leaving the numerators in (Ot-, PT^j alone, 
and discretizing the denominators: 

•k t Xtix, 


On = ~- 


Qn(X) 


24 











and 


l rT, 


Pn(£) = -- 


fo n Vt n W? 


Qn(X ) + (^) 2 Q„(S) 

where 

= X t + 9 n X t , 0 ^t^T n . (45) 

and we recall that Q n (Z) is a notation for the Riemann-sum rectangle approximation ^ ) 2 ■ 

We also consider the version of p n (£) based only on discrete observations of X: 


£ /o^ Wt" 


Pn(£) = ~ 


Q„(X) + (6» n ) 2 Q n (S) 


Combining Lemma [3] and the almost-sure convergence (j37|) we deduce the strong consis¬ 
tency of Q n . 

Theorem 12 Assume H £ (1/2,1). If A n ^ n a /or some a € (—oo, 1/H), then 

6n —> e* 


almost surely as n —>• oo. 

By Lemmas l6l and fill it is easy also to deduce the strong consistency of p n (S) and p n (X). 
Theorem 13 Assume H € (1/2,1). If A n ^ n a for some a € (—oo, 1/H), then 

PnV) —>• P* 


almost surely as n —» oo. 

In addition, if n 1+a A% +1 —>• 0 /or some a > 0, 

p n (S) —► P* 


almost surely as n —>• oo. 

To establish the asymptotic normality of (0 n , p n 1, we can write 







e*y/T n {^ - Q n (X)) 
Qn(X) 


Similarly, 


y5(?» - pi = - /) + AlAXXiSA , 


Theorem [8] provides the convergence of the last summands in each of the two lines above. 
Combining this with the convergences we obtained in Lemmas [12] and [IT] we obtain the 
following result. 
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Theorem 14 Let H € ( 5 , |) and nA ^ +1 —> 0. Then 

(£„ - e*,p n {T) -p*)^Af (0,* p r P) 

In addition if n 3 A 2 n H+3 —>• 0, 

(X - 0*, Pn&) -p*y^N (0,* P T P) 
where P the matrix defined in < 1351) . 


4.2 X and E are observed 


The problem with the auxiliary estimators Q n and p n is that they still contain Skorohod 
integrals. In order to devise a further scheme that allows us to evaluate them, at least ap¬ 
proximately, using discrete data only, we begin by using the discrete observations of X and 
E, and recalling that, from Lemmas [3] and [6] we have 

Jo Xt6Xt —* -p°(p + e )v x , 

J Jo ^ t5Xt —* -P 0 (P + 9 )P^ 

where r) X and ?/ s , which are also functions of H.O^p, are given in Lemma [2j Since these 
limits depend on the parameters we are trying to estimate, one strategy is to rewrite the 
strong consistency results of Theorems fl2l and 1131 for Q n and p n (E) as implicit definitions of 
new estimators, where the numerators in the definitions of 9 n and p n (E) are replaced by their 
limits recalled above, and each instance of 8 and p therein are replaced by the new estimator 
we are trying to define. The same substitution must be done with the expressions 8* and p*, 
since these are the limits of 8 n and p n . In other words we consider only that the denominators 
in 9 n and p n contain data, and replace all other instances of ( 8 , p) in the limits in Theorems 
m and [13] by the pair of estimators we are trying to define. After some minor manipulations, 
this leads to the following definition of a new pair of estimators (0 n ,/5 n ) as solution of the 
system of the following two equations, if it exists: 

h , x _ 

°n + Pn- (p n -0 n )Q n (X) 

( _ ( Sn) 2 -(pn ) 2 _ _ _ HT(2H) _ 

[( 9r l ) 2 - 2H -(pn) 2 - 2H + {p n +en)\(p n )- 2H -(e rl )-MT)J _ Qn (X) + (8„ +Pn ) 2 Qn (E) 

We emphasize that the above is an implicit definition of (0 n ,p n ). It is also rather opaque. 
The system (1461) can be simplified slightly using more elementary manipulations. We find that 
the definition of ifimPn) is equivalent to the following: 



F(e n ,p n ) = (Qn(X),Q n ( E)) 

where F is a positive function of the variables (x, y ) in (0, + 00) 2 defined by: for every (x, y ) € 

( 0 , +00) 2 


F(x,y) = HT(2H) x 


^ ( y 2-2 H _ x 2-2H^ x -2H _ y -2H) if x + y 

((1 — H)x ~ 2H , Hx~ 2H ~ 2 ) if x = y. 


(47) 
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Interestingly, this shows that a good candidate for the discrete version of the least-squares 
estimator of (6, p ) is none other than a type of generalized method of moments estimator 
obtained via Lemma[2]after discretizing the expressions St(X) : = TX 2 ds and Sp(T,) : = 
T~ l Y? s ds. We now consider the question whether System (l46l) has a unique solution 
(# 2 i/5n)> and how this may imply strong consistency for these estimators. 

Since for every ( x , y) E (0, +oo) 2 with x y 


J F ( x,y) = T(2H+1) 


'(1 -H)x 1 - 2H (x 2 -y 2 )-x(x 2 - 2H ~y 2 ~ 2H ) 

(x 2 -y 2 j 2 

Hx- 2H - 1 (x 2 -y 2 )+x(x- 2H -y ~ 2H ) 
(x 2 -y 2 f 


(1 -H)yl- 2H ( y 2 -x 2 )-y( y 2 - 2 X-x 2 - 2 *) 
(x 2 -y 2 ) 2 

Hy- 2H ~ 1 (y 2 -x 2 )+y(y- 2H -x~ 2H ) 

(x 2 -y 2 ) 2 


the determinant of Jp ( x , y) is non-zero on in (0, +oo) 2 . So, I 7 is a diffeomorphism in (0, +oo) 2 
and its inverse G has a Jacobian 


J G (a, 6) 


T(2H + 1) 
det Jp (x, y) 


Hy- 2H - 1 (y 2 -x 2 )+y(y- 2H -x- 2H ) (1 (y 2 -x 2 )-y{y 2 ~ 2H -x 2 ~ 2H ) 

( x 2 — (x 2 — 

Hx~ 2H ~ 1 (x 2 —y 2 ^+x(x- 2H —y- 2H ^) ( 1 —H)x 1 ~ 2H (^x 2 —y 2 ^—x{x 2 ~ 2H —y 2 ~ 2H ^j 

( J 2 -y 2 ) 2 (x 2 -y 2 ) 2 


where ( x,y ) = G(a,b). 

Hence, (1571) and ()42[) lead to 

{On,Pn) =G(Q n (X),Q n (E)) — >G(?] X ,p E ) = (e,p ) 


almost surely as n —>• oo as soon as A n ^ n a for some a E (— 00 ,1 /H). Summarizing, we have 
proved the following. 


Theorem 15 Let H E (1/2,1) and assume that A n ^ n a for some a € (— 00 ,1 /H). Then, 
as n —> 00 

( finiPn ) t ( 9,p ) 

almost surely. 

We may now prove a normal convergence result for (d n , p n ) based on Theorem fill Note 
that the second part of Theorem Q31 is not needed here because we rely on fully observed E in 
this section. 


Theorem 16 Suppose that H E (I, |) and nA^ +1 —> 0. Then 

VT~n{e n -e,p n -p) A/"(0/ M f P T P M) 

where the matrices T, P and M are defined respectively in HA-5\) . \35 1) and 11501) . 
Proof. We have 

0 n = 

and 


k Io n X t SX t __ J0 {X) 


Qn{X) 


Qn(X) 


x JX yn 5 yn 


Pn(E) = ~ 


4(H) 


Q n (X) + (0 n ) 2 Q n (Z) ’ Q n (X) + (0 n ) 2 Q n (E) 
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Then, we can write 


\fr~n (o n — o, p n — p)) 


= Vr~n G 


( rl ( J n( x ) ,2^-2 (j^V) J 9 n (X)' 




On 


( 0 : 


yPn(^) 0 n 


- ( o,p ) 


= \fr~n (GoL (j 9 (X), JP(V),O n ,p n (T,)) - (0,p 


where L(r,s,u,v) = (-, -4 - 4)- 

On the other hand for any e € (0,1) 


T rt 


fl 


r 2H ~ 2 e~ r drdt = T(2 H - 1) + o(—) 

T £ 


because 


1 t'T roo i /*oo 

1 / / r 2H-2 e -r drdt ^ ^_ { l_ e -T/2 ) ^H-2^/2^ 

JO 


J'l—e 


ff 


2T l ~ £ 

0. 


Combining this together with (I23[) and the point 5) of Lemma [21] we can write 


1 


Jn(X) = (p + 0)p X + o{—=) 

V J-n 

Lm> 


(48) 


where o(-^=) denotes a random variable such that fTlo(-^=) converges to zero almost surely 
as T n —>• oo. 

Similar argument leads to 


1 


Since 


we can write 


jp(y) = pe{p + 0)^ + o{-=). 

v -L n 

GoL((p + 0)r, x ,p0(p + 0)r,*,6*,p*) = ( 6,p ) 


(49) 


V^n (On — 0,p n — p)) 


= sjT n (GoL ( J° n (X), JP(V),O n , Pn(S)J - GoL ((p + 0)r, x ,p0(p + 0) V E ,0*,p 

= ^[GoL (j e n (X),J£(V),0 n ,Pn(Z)) - GoL ^p + 0)p x ,p0(p+ 0)^,0^^)) 

((P + 0)v x , —p0(p + Off, 0 n , Pnfjfj - GoL ((p + 0)p x , p0{p + Q)rf, 0*,p *) 


+ GoL 

Sn T T'n • 


From (1481) and (l49l) we obtain s n —)• 0 almost surely as n —>■ oo. 
On the other hand, by Taylor’s formula 


Tn = VTn 0 n - 0*,p n ( E) - p* ) M + d, 
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where 


M = 


(50) 


with 




!?(«' 

§?(«■ 


* 

>P 


h(u,v) = (hi,h 2 )(u,v) = GoL ((p + 6)i] X , pd(p + 6)^,u,v) . 


We can write 


h(u,v) = (Gi,G 2 )og(u,v) 


where 

= { 9 i,g 2 )(u,v) = L ((p + 0)i] x , p6(p + 0)r]-, u,v) . 

Moreover for i = 1,2 


/ \ 9Gi , , \ \ / \ . &Gi / / \\ ^72 / \ 


and 


dhi dGi ^ d 9i, \ , dG i ( r \\ d 92, , 

^(u,v) = ^( S (u, v ))-^(u,v) + ^( 9 (u,v) ) — (u, v) . 


On the other hand, d n converges in distribution to zero, because 


IKII < c{H,6,p)Vr^\\(Pn(x) - 0*,e n - P *) II 2 . 


It is elementary that if for any u € Cl there exists no(w) € N such that X n (uS) = Y n (u ) for all 

n ^ no(w) and X n Oasn^ oo, then Y n 0 as n —>• oo. 

Combining this with Theorem 1141 the proof is completed. ■ 


4.3 X is observed 

In the previous section, we encountered theorems in which X and £ are both assumed to be 
fully observed in discrete time. Since £ is the time-antiderivative of X , such an assumption 
corresponds, for instance, to the physical situation where X is the velocity of a particle, and 
£ is its position. 

In this section, we abandon such a framework, and assume instead that only X is observed 
in discrete time. Thus we consider the following pair of estimators (j) n ,PnJ : 

( e n ,Pn ) =G(Q n (X),Q n (£)) 

where the deterministic explicit function G was identified in the previous section as the inverse 
of the function F given in (1471) . Equivalently, is the solution of the system (|46l) . or 

its equivalent form (1471) . with £ replaced by the process £, which relies only on observations 
of X. Using same arguments as in Section 14. 2 1 and Lemma HT1 but relying now on the second 
part of Theorem [13] (hence the stronger condition on A n for the strong consistency result) 
and the second part of Theorem Q31 (hence the stronger condition on A n for the convergence 
in law result), we conclude the following. 
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Theorem 17 If n 1+a A^ +1 — > 0 for some a > 0, 

{Oni Pn'j t {91 P) 

almost surely as n —>• oo. 

Theorem 18 Let H € (^, |). If n 3 A^ H+3 —> 0, then, as n —>• oo 

where the matrices T 7 P and M are defined respectively in HA-5\) , i t 35)) and K50\) . 


4.4 X and U are observed 


When both X and V are observed, the estimator of 9* based on continuous data is 6t given 
in Q but the estimator of p becomes the usual full-observation estimator of an fBm-driven 
Ornstein-Uhlenbeck process as in M, i-e. 


Pt ~ 


Io T vm 


lo V?dt 


Following similar arguments as in the beginning of Section 14.21 the natural candidate for the 
estimator based on discrete data of X and V is the pair (jl n , p^j defined as the solution of 
the following system: 


+ P n ~ 

( HT(2H)\ 2 h 
P n ~ V Qn(V) ) 


HT(2H) r ,_ \2—2H _ (n \2-2H] 


Qn{X) 


We see that p n is defined explicitly autonomously via the discrete-data-based statistic Q n (V). 
With p n now known, elementary manipulations yield that 9 n is precisely the solution of the 
following simple equation 


(Or 


\2-2 H ( Qn(X ) 


\HT(2H) 


(oj 2 = (p n ) 


2~2H / Q n ( X ) \ 


\HT(2H)) \— 1 


P. 


Define 


y2- 2H - x 2- 2H 2 H 

y 2 —x 2 ’ & 


if x ^ y 
((1 — H)x~ 2H , if x = y. 


F(x, y) = HT(2H) x 

l U 1 - H 

Its Jacobian is given, for every (x, y) € (0, +oo ) 2 such that x 7 ^ y, by 
Jp(x,y) = r(2P+i) 


(1— H)x 1 2H (x 2 —y 2 )—x{x 2 2H — y 2 2H ^j ( 1 —H)y 1 2H {y 2 —z 2 )—y(y 2 
(x 2 -y 2 ) 2 (x 2 -y 2 )' 2 


0 


Thus the Jacobian of is inverse G is as follows 


-y 


-2H-1 


2 - 2 H ~2 — 2H \ 


Jq ( a,b ) — 


T(2H + 1) 
det J T (x, y) 


-2H-1 (1 -Hjy 1 2H (y 2 -x 2 )-y(y 

(1 -H)x^~ 2h {x 2 -y 2 )-x{x 2 - 2H -y 2 ~ 2H ) I ’ ^ = G ^ ' 


{x 2 -y 2 ) 2 


Using same arguments as in Section 14.21 we obtain 
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Theorem 19 Assume that A n ^ n a for some a G (— 00 , 1 /H). Then, as n 


00 


{On,P n ) -> (°iP) 

almost surely. 

Theorem 20 Suppose that H G (^, |) and nA^ +1 —>• 0. Then, as n —>• 00 


(«n -6,p n -p)^ AT^Of Q l PT PQ) 


where Y and P are defined respectively in HA-A) and USA) , and where 

'§(£*,p) ^ 


Q = 


K o*, P ) 
(o\ P ) 


such that 

with 

and 


(fi,f 2 )(u,v) = Gol (u,v) 

l{u,v) = L((p + 9)i] X ,p 1 - 2H ,u,v) 

L(r,s,u,v) : = () . 

\u vJ 


5 Appendix 

In this appendix, we present some calculations used in the paper. 

Fix T > 0. Let f,g : [0,T] — > K be Holder continuous functions of orders a G (0,1) and 
(5 € (0,1) respectively with a + ft > 1. Young |29] proved that the Riemann-Stieltjes integral 
(so-called Young integral) f 0 f s dg s exists. Moreover, if a = /3 € (^,1) and (j) : R 2 —> R is 
a function of class C 1 , the integrals f 0 ^y(/ u , 9u)df u and f 0 9u)dg u exist in the Young 

sense and the following chain rule holds: 

ft deb ft deb 

<t>(ft,9t) = J>(fo,9o) + / (fu,9u)df u + -w-(fu,9u)dg u , O^t^T. (A-l) 

Jo of Jo og 

As a consequence, if H G (^, 1) and [ut, t G [0,T]) is a process with Holder paths of order 
a G (1 — H, 1), the integral u s dB^ is well-defined as a Young integral. Suppose moreover 
that for any t G [0, T], ut G D 1 ’ 2 , and 

P (^j J \D s ut\\t — s\ 2H ~ 2 dsdt < oo^j = 1. 

Then, by [1], u G Domd and for every t G [0, T], 

Us dB* = [ u s 6B” + H(2H - 1) [ [ D s u r \s - r\ 2H ~ 2 drds. (A-2) 

Jo Jo Jo 
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In particular, when tp is a non-random Holder continuous function of order a € (1 — H, 1), we 
obtain 






(A-3) 


In addition, for all ip, if> € \H\, 

( ! I I ! I I . ! I I ' J I 

f p s dBj? f ipsdB^i = H(2H — 1) f f p{u)if(v)\u — v\ 2H ~ 2 dudv. 

Jo Jo J Jo Jo 

Lemma 21 Let m, ml > 0 and let X m be the process defined in 1151) . Then, 

1) A (m, m') := H(2H— 1) / 0 °° e~ ms e~ m ' r \s—r\ 2H ~~ 2 drds = + m' l2H ^j, 

2) 0 ^ A(m, m') - E{X™X™') < c{H, m, m')e-*/ 2 , 

V sup t> 0 £[|X t m | p ] ^ c(H, m,p) < oo, 

4) 0 ^ ^ c(H,m,m')\t- s\ 2H ~ 2 , 

X 171 

5) For every e > 0, ^—> 0 almost surely as T —>• oo, 

6) E{\X™ - X™\p) ^c{H,m,p)\t- s\p h . 

Proof. To prove equality 1), we just write 


poo poo 

X(m, m!) = H{2H-1) / / e~ ms e~ m ' r \s - r\ 2H ~ 2 drds 

J 0 J 0 

poo ps 

= H(2H-1) dse~ ms dre~ m ' r \s - r\ 2H ~ 2 

Jo Jo 

poo poo 

+H{2H-1) dse~ ms / dre- m ' r \s - r\ 2H ~ 2 

Jo Js 

HT(2H) f j —2H , ri-2H\ 

= - - [m + m 

m + m V / 

For the point 2) see [15] . For 3) and 6) we refer to fl6] . and for 4) and 5) see [14j Lemma 5.2 
and Lemma 5.4] ■ 

Theorem 22 Let H € ( 5 , |). Define for m > 0 

Then, as T —> 00 , 


-j=(h(fr)Mft)) ^AA( 0 ,r) 


(A-4) 


where T is a symmetric nonnegative definite matrix which has the following explicit expression 


r = r)H 


h h 
h h 


(A-5) 
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where h = 6 1 4H , l 2 = p 1 4H , h = (ah-i^-O'A 4H ~ P 1 4/i ] and 


m = H\AH - 1) 


T(2H) 2 + 


T(2H)T(3 - AH)T(4H - 1 ) 

r(2 - 2H) 


Proof. Notice that for (1A-4|) to hold it suffices that prove that for every a, b £ M, 


G T :=aI 2 (ft) + bI 2 (f p ) 


converges in law to M (0, (a, b)T t (a, b )) as T —> oo. 

Fix a, b € R. Since Gt is a multiple integral, by the isometry property of double stochastic 
integral I 2 , we get the variance of Gt as follows 


EGt = 


a 


H 


ri 

,2 1 T 


T 2 T 3 

+ b 2 !z + 2ab^ 

' rj~\ ' nrt 


where 


4 = 


7 ' 2 - 
1 T — 


e -6\t-s\ e -e\u-v\\ t _ n |2^-2| s _ v fH-2 dtdsdudVj 

'[0 ,T] 4 

[ e -p\ t - s \ e -p\ u ~ v \ 1 1 - u\ 2H - 2 \s - u| 2 H " 2 dtds(iu(ir;, 
«/ [0,T] 4 


To _ 

irp - 


_ d 2 ^- 2 |s - d 2 ff - 2 didsdlKfo. 


' [0,T] 4 


Using the same argument as in the proof of III Theorem 3.4], we have as T —>• oo 

= mh 


a 2 H 4 
lirriT- H T 


and 


2 T 


a 2 H 4 

lirriT-too = r lH^2- 


Now, let us estimate I? r . We have 


Hr 


= 2 


' [0,T] 3 


e ~p(T~s) e -6\u-v\ (ji _ u )2H-2| s _ ^|2H-2 dsdud?; 


[ e -p\ts\ e -<>(T-v)( T _ t )2fT-2| s _ v |2^-2 dsdudu 

j ro.Ti 3 


+2 

'[0,T] 

:= 0 ) + A T ($, p). 

Making the change of variables T — s = x , T — u = y and T — v = z 

A T (p, 6) = 2 [ e- px e- ely - zl y 2H ~ 2 \z - x\ 2H - 2 dxdydz. 

J [0,T] 3 

This implies 

UmT^rooAT(p, 0) = Aoo(p, 9) = 2 f e~ px e~ 9 ^ z ^y 2H ~ 2 \z — x\ 2H ~ 2 dxdydz. 

J [0,oo] 3 
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Making the change of variables z — x = w, we obtain 


4x> (p, 0 ) = 2 


= 2 


*[0,oo) 2 J—x 
r ry-x 

' [0,oo) 2 J—x 


+2 


/• /*oo 

«/ [0,oo) 2 J —X 

e - pXe -e(y-w-x) y 2H-2\ w \2H-2 dwdxdy 

P POO 

/ / e~ px e d ^ y ~ w ~ x ' > y 2H ~ 2 \w\ 2H ~ 2 dwdxdy 

J ('O.oo') 2 J v—x 


Integrating in x we get 


Aoo{p,9) — 


(-p + 0 ) L 


/* OO /»-)-00 

/ / ( e (-P+e)( 2 /-^) _ e (-P+ 0 )[(-^) v °])l [(2/ _ j/;) _ (( _ 

VO J— OO 


id)V0)] 4 


xe w )y 2H 2 \w\ 2H 2 dwdy 


{-P-0) 
2 


oo />+oo 


n 


o(-0-p)[{y-w)V 0] 0(y—w) 2H—2\ \2H—2 


^y-wiyZM-z\ w \^^dwdy 


-A 1 (p, 8) H ——A 2 (p, 9). 
9 -p p + 9 


Furthermore 

A 


POO POO 

\p,8) = / (e~ p ( y+w) -e- pw e- 6y )y 2H - 2 w 2H - 2 dwdy 

J 0 J 0 

/•oo /•?/ 

+ / / {e-* y - w) - e- e{y - w) )y 2H - 2 w 2H - 2 dwdy 

Jo Jo 

( /•oo \ 2 / /-oo \ / /-c< 

L yW ~ 2e ~ ndv ) - (i y 2H ~ 2 e ~ yy djj (J" 

POO POO 

+ / [ e -p(y~ w ) - e - e ^- w) )y 2H - 2 w 2H - 2 dydw 

J 0 J w 


w 2 H-2 e -pw dw 


' 0 J w 

= r (2 H - lfp 1 - 211 [p 1 - 211 - 9 1 - 211 ] 

POO POO 

+ / (e _pa: - e- fe )(x + w) 2H - 2 w 2H ~ 2 dxdw. 

Jo Jo 

Using (x + w) 2H ~ 2 = r ( 9 ^ 2 g) fo° Z 1 ^ 211 e~^ w+x ^d£, the term A 1 becomes 

/ 4 1 (/o, 0) = r(2 h - 1 ) 2 P l ~ 2H [p l ~ 2H - e 1 - 211 } 

T(2H-1) ro ° ro ° 


+ 


r(2 - 2 H) 


+- 


POO POO 

/ / Z 2 - 4H e-t x {e- px - e~ 9x ){d^dx 

Jo Jo 

[P 2 


= v(2H -1 f P l ~ 2H [p l ~ 2H - e 1 - 211 ] 

T( 2 H - l)r (3 - 4H)T(4H - 2) , 2 _ 4H q2 _ 4H 


r(2 - 2 H) 
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Similarly, we obtain 


pOO f'OO 

A 2 (p, 9) = e-^ v+w) y 2H - 2 w 2H - 2 dwdy 

J 0 Jo 

poo py 

+ / e- p ( y - w ^y 2H - 2 w 2H - 2 dwdy 

Jo Jo 

poo poo 

+ / / e 9 ^y 2H - 2 w 2H - 2 dwdy 

Jo Jy 

= p 2 ~ 4H T(2H - l) 2 

T(2H - l)r(3 - AH)T(AH - 2) ^ 2 _ 4 „ + q2 _ ah ^ 


+ 


Thus, Aooidjp) is also obtained. 
Consequently 


T (2 - 2 H) 


II 


UrriT^ooatHY = a H( A oo(p,0) + A oo (0,p)) 

= 2 ijiih- 

Finally, combining the above convegences we deduce that as T —> oo 

EG 2 t —> (a, b) T *(a,6). 

On the other hand 

D s G t = -j= (J* (ae~ 6 ^ + be-^ 3 -^ 5B t + £ (ae~ e ^- s) + be~ p{t ~ 

a,b I \r a ib 


Vt 




Hence 


s,T ' 


\\DG t \\h = J Q J Q 

■ £ £ (x? b x? b + 2X?’ b Y s a * + \r-s\ 2H 


X a s ' b + Y ^ b J (X r a ’ 6 + Y r a *J |r - s| 2 ^" 2 dsdr 


T 


a.,6 , D a,& i /~ia,b 

—— ( Arp + + Crp, 


Since Xg’ b belongs to the first Wiener chaos of B H 
2 


E 


A^ b - EA^ b 


= 2 E(X“’ b X?’ b )E(X“’ b XZ’ b )\u - r\ 

2[0,T] 4 


2H ~ 2 \v - so¬ 


using similar arguments as in III Lemma 5.4 of web-only Appendix], 


A^ b - EA^ b 


€ c 


( h,o,p ) [ 

J\ c 
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r p4—SH 


[0 ,T] 4 


% 1] 4 


I \2H—2 1 1 2H—2 1 \2H—2 1 
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s — r u — it u — r \v 



2 dsdr 

2 dsdrdudv. 

- s\ 2H ~ 2 dsdrdudv 

- s\ 2H ~ 2 dsdrdudv. 
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Using the same argument for B^’ b and C^ b we conclude that 




E {\\\DG t \\h - E\\DG t \\ 2 h 
c(H,d,p) r 


2 I 2 


2^6-8 H 


\2H—2 1 1 2H—2 1 \2H—2\ \2H-2 > r , r 

s — r\ \v — u\ \u — r\ \v — s dsdrdudv 


'm 4 


0 


as T —> oo, because H < This completes the proof of Theorem 
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